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THE EFFECT OF WALL INTERFERENCE UPON THE AERODYNAMIC 
CHARACTERISTICS OF AN AIRFOIL SPANNING A 
CLOSED-THROAT CIRCULAR WIND TUNNEL 
By Valter 0. Vincent 1 and Donald J. Graham 

SUMMARY 


The results of a theoretical and experimental Investi- 
gation of wall Interference for an airfoil spanning a closed- 
throat circular wind tunnel are presented. Analytical equa- 
tions are derived which relate the characteristics of an air- 
foil in the tunnel at subsonic speeds with the characteristics 
in free air. The analysis takes into consideration the ef- 
fect of fluid compressibility and is based upon the assumption 
that the chord of the airfoil is small as compared with the 
diameter of the tunnel. The development is restricted to an 
untwisted, cons tant— chord airfoil spanning the middle of the 
tunnel. Brief theoretical consideration is also given to the 
problem of choking at high speeds. Results are then presented 
of tests to determine the low— speed characteristics of an NACA 
*4412 airfoil for two chord— diameter ratios. While, on the 
basis of these experiments, no appraisal 1b possible of the 
accuracy of the corrections at high speeds, the data Indicate 
that at low Mach numbers the analytical results are valid, 
even for relatively large values of the chord— d iamet er ratio. 


INTRODUCTION 


The design of modern high-performance airplanes requires, 
insofar as possible, an accurate knowledge of airfoil profile 
data at Reynolds and Mach. .numbers attained in flight. Since 
the size and power of wind tunnels are subject to various 
practical limitations, most existing tunnels, even if they can 
provide the desired Mach number, are not capable of attaining 
full-scale Reynolds numbers for all flight conditions. To 
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minimize this shortcoming in tunnel tests of airfoil profiles, 
it is therefore necessary to use models having as large a 
chord as possible relative to the cross-sectional dimensions 
of the tunnel test section. In order to eliminate the effects 
of supporting struts and to exclude the indeterminate tunnel- 
boundary interference involved in the testing of large— chord 
airfoils of limited span, it has become common practice in 
euch tests to use airfoils which completely span the test sec- 
tion. Even for these so-called "through" models, however, the 
tunnel— boundary interference can still be considerable, and 
accurate correction must bo made for its effects if the tunnel 
data are to be usod with confidence in tho calculation of free- 
fligiit airplane characteristics. 

The tunn el— boundary interference for airfoils spanning 
wind tunnels of various types has been the subject of numerous 
theoretical and experimental investigations. The Interference 
for rectangular tunnelB having rigid walls normal to the span 
of the airfoil and either rigid walls or free boundaries paral- 
lel to the span has been discussed theoretically by several 
writers. Tor example, Lock (reference l), Glauert (reference 
2), and Goldstein (reference 3), give the necessary tunnel— wall 
corrections for an airfoil spanning a rectangular tunnel in an 
incompressible fluid; ^hilo Goldstein and Young (reference 4) 
show how these corrections, as well as those for any general 
case of interference m an incompressible fluid, can be modi- 
fied to take account of fluid compressibility. Eeference 5 
gives the corrections for the compressible case in a closed- 
throat rectangular tunnel, as well as a critical discussion 
of the results of the previous references and some experi- 
mental data from low— speed tests. Tage (reference 6) also 
presents experimental drag data for several symmetrical bodies 
of various sizes In a clos ed— throat rectangular tunnel. Ex- 
perimental and theoretical results for an airfoil spanning a 
completely open— throat rectangular tunnel are given by Stiiper 
(references 7 and 8). The case of an airfoil spanning an 
open— throat circular tunnel has been the subject of a number 
of investigations, including theoretical treatments by Glauert 
(reference 9), Stuper (references 7 and 8), and Souire (ref- 
erence 10), and experimental measurements by Stupor (refer- 
ences 7 and 8) and Adamson (reference ll). Apparently, the 
case of the closed— throat circular tunnel has received no at— 
t ent ion . 

Since this last case is ofton encountered in practice-, 
an investigation was made of the tunnel— wall interference at 
subsonic speeds for a wing spanning a closed— thr oat circular 
tunnel. The present paper presents the results of this in— 
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vestlgatlon. In the first part of the paper, analytical equa- 
tions are derived relating the characteristics of the airfoil 
in the tunnel with those" In free airfor a compressible fluid. 
Some consideration is also given to the phenomenon of choking 
which occurs at high speeds. In the second part, the validity 
of the theoretical results is examined by the analysis of ex- 
perimental data for an HACA 4413 airfoil f»r two ratios of 
airfoil chord to tunnel diameter. The investigation is re- 
stricted to untwisted constant— chord airfoils spanning the 
middle of the tunnel. 


THEORY 


As in reference 5, the theoretical development of the 
tunnel— wall corrections is divided conveniently into two gen- 
eral sections. First, the Influence of the wall upon the 
field of flow at the airfoil in the tunnel is determined. 
Second, the aerodynamic characteristics of the airfoil in this 
field of flow are related to the corresponding quantities in 
free air. In thiB manner, simple formulas are finally ob- 
tained which enable the prediction of the free— air character- 
istics when the characteristics in the tunnel are known. 

Again as in reference 5, the analysis is based upon the 
method of superposition. To this end, it is assumed that the 
airfoil is of small thickness and camber, so that the Induced 
velocity is everywhere small as compared with the velocity, 
of the undisturbed stream. With this assumption, the total 
Induced velocity at any point is the simple vector sum of the 
separate velocities Induced at that point by the interference 
between the tunnel wall and the airfoil camber, thickness, 
and wake. Thus the effects of camber, thickness, and wake 
may each be analysed separately and superposed to obtain the 
desired result for the complete airfoil. Ab pointed out in 
reference 5, this procedure is permissible even in the com- 
pressible fluid if the airfoil is of small thickness and cam- 
ber as assumed. 

Before proceeding to the actual development of the theory, 
it is useful to contrast the present problem with the problems 
of through airfoils in the various types of rectangular tunnels 
and in the open— throat circular 'tuhnel^' ■ In^^lre -eas-eof an 
airfoil spanning a rectangular tunnel having rigid wallB nor- 
mal to the span of the airfoil, the problem is relatively 
simple. If the effect of the boundary layer along the tunnel 
walls is neglected, the flow is sensibly the same in all planes 
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normal to the span; that la, there 1 b clearly no spanwise 
.variation in lift. The air flow is thus essentially two- 
dimensional. and the Interference nrohlems of camber, thick- 
ness, and wake can he analysed by the customary means of a 
system of images with axes parallel to the span of the air- 
foil (references 1, 2, 3 » an< ^ 5 )* This is true whether the 
tunnel boundaries -oarallel to the span are fixed or free. 

In this manner, tunnel-boundary corrections can be derived 
for airfoils of moderately large chord as compared with the 
height of the tunnel test section. 

In the case of an airfoil spanning a completely free 
jet, whether rectangular or circular in section, the lift 
necessarily falls to zero at the boundary of the jet. There 
thuB exists in this case a pronounced spanwise variation in 
lift and an attendant system of trailing vortices. In the 
existing treatments of the problem, only the interference 
between these trailing vortices and the jet boundaries is 
considered, the interference effects associated with the 
ehordwise distribution of bound vortices and with the air- 
foil thickness and wake being completely neglected. This 
procedure implies the assumption that the chord of the &ir~ 
foil is very Bmall relative to the dimensions of the jet. 

In this manner, the problem is reduced to a limiting case 
of the usual problem of an airfoil partially spanning the 
jet, and, as in this latter case, the component of down- 
wash induced at the airfoil by the interference between the 
walls and the trailing vortices is one-half as great as the 
corr es-nonding component an infinite distance downstream. 

The theoretical determination of the wall interference may 
thus be treated as a problem of two-dimensional flow in a 
plane normal to the axis of the tunnel infinitely far behind 
the airfoil. The boundary conditions for either the 
rectangular or circular jet are then readily satisfied by 
the introduction of a suitable system of image vortices with 
axes parallel to the axis of the tunnel (references 7» 8, 

9, and 10). This method of analysis, however, is inadequate 
if the chord of the airfoil is even moderately large as 
compared with the dimensions of the jet. 

The case of the airfoil spanning a clos 6d-throat circu- 
lar tunnel is more complex than either of the foregoing prob- 
lems. Unlike the condition prevailing in the free jet, the 
lift in thiB case need not fall to zero at the boundary - 
that is, at the tunnel wall - so that the spanwise variation 
in lift is not necessarily large. In fact, aa will be seen, 
the lift is constant across the span of the airfoil, and no 
system of trailing vortices exists. The assumption of a very 
small chord and the consequent reduction of the -problem to a 
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case of two— dimenaibnal flow In a plane infinitely far down- 
stream i® thus without meaning. On the other hand, an analy- 
sis for airfollB of moderately large oh'ord ' in the manner em- 
ployed in the case of the rectangular tunnel with rigid side 
walls is not possible. In the clos ed— throat circular tunnel 
the flow in all planes normal to the span of the airfoil is 
not the same, so that the effect of the hound vortices, and 
of the airfoil thickness and wake as well, cannot he treated 
as a problem in two-dimensional fl&w. furthermore, the bound- 
ary conditions at the tunnel wall cannot be satisfied for the 
actual three-dimensional problem by any known system of Images. 
The solution of the problem for the closedr-thr oat circular 
tunnel thus requires an analysis entirely different from those 
employed in the previous instances. 


Influence of Tunnel Vail upon field of flow at Airfoil 

An approach to the problem of the airfoil spanning a 
closed— throat circular tunnel is afforded by the work of Ton 
Xarman and Burgers in reference 12 ( pp . 266 to 273), where 
the velocity potential at an arbitrary point in the tunnel is 
determined for a U— shape vortex of infinitesimal span in an 
incomer ess ible fluid. 

A system of rectangular coordinates- x, y, z 'is intro- 
duced as shovn in figure 1. The x— axis is taken on the center 
line of the tunnel with its positive direction downstream. 

The z— axis is positive downward, and the y^-axis positive to 
the left for an observer looking against the direction of flow. 
An alternative Bystem of cylindrical coordinates x, u> , 0 is 
defined by the relations 


y 

z 


U) 008 B 

u; sin 0 J 


( 1 ) 


The positive direction of circulation is defined so that a 
vortex with positive circulation exerts a ferce on the fluid 
in the direction of the positive a— axis. In other words, 
the lift force experienced by a positive vortex is in the 
pegative z direction. The velocity of the fluid in the un- 
disturbed stream is denot ed~ by J V 1 ' And the radius of the 
tunnel by r. Other symbols are defined as introduced in the 
text. A list of the more important symbols and their defini- 
tions is given in appendix C. 
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Consider now a TJ-fshapa vortex of infinitesimal span dp 
parallel to the y>-axia and situated in the yz— plane at tho 
point p = u> 0 cos 0 O , $ = u> 0 sin 0 O . If the strength of the 

vortex is denoted hy I' * , the velocity potential in the 
closed tunnel at the points x, U) , 0 is given by Ton Karman 
and Burgers, for negative values of x, as 


OB 



( Jt should be noted that the quantity I appearing in these 
equations is merely a variable of integration and has no 
physical significance.) The quantity J n (\ a u:) is a Bessel 
function of the first kind of the order m. The summation 
with respect to m extends over all the positive integers 
and includes m = 0; the prime added to the summation sign 
indicates that a factor l/2 must be inserted before the 
term corresponding to a = 0. The summation with respect to 
8 for every m extends over all positive roots of tho equa- 
tion 

J n ? (V) = 0 (4) 

where .Jn (X. s r) ie the derivative of the function J a (X. 8 r) 

with respect to its argument. The notation used throughout 
this paper for tho Bessel functions is that of Watson (refer- 
ence 13;, which is the same as that of the Smithsonian tables 
(refer ence 14 ) . 

By differentiating ft with respect to £ and then in- 
tegrating with respect to £ as indicated in equation (s) f 
the velocity potential becomes finally 
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r»dn 


00 oa 

1 X 

m= 0 b c 


e XflX J n (A B u>) 


ri . C 1 -^) x » aj " au ' r) 

X f"m cob 0o sin m( 6—6 0 ) — B — 0 - 

L ®o ■ 

+ X. s ain0 o cob mCB— 0 O ) Jj/CXbUJo) J 


(5) 


As pslntsd out, this expression applies only at negative val- 
ues of x. As will be seen later, the necessary results for 
positive values of x can he derived from considerations of 
symmetry . 


.By means of equation (&), it is possible to evaluate the 
wall Interference associated with both airfoil camber and 
thickness for the case of the incompressible fluid. These 
results can then be modified for the effect of fluid compres- 
sibility by the methods of reference 4. It is found finally 
that, for a closed— throat circular tunnel, the effects of in- 
terference between the walls and the airfoil camber are iden- 
tical with the corresponding interference effects for the 
same airfoil spanning a closed— throat rectangular tunnol , the 
height of which bears a known relation to the diameter of the 
circular tunnel. A similar conclusion is obtained regarding 
the effects of interference between the vails and the airfoil 
thickness, except for a numerical difference in the relation 
between the diameter of the given circular tunnel and the 
height of the equivalent rectangular tunnel. The interfer- 
ence effects associated with the wake of the airfoil are not 
analyzed in detail, but their magnitude can be estimated with 
reasonable accuracy by comparison with the results for the 
thickness effect. In order to simplify the complex mathe- 
matics of the problem, the Interference effects are calcu- 
lated only for the section of the airfoil at the center line 
of the tunnel. As will be seen later, however, experimental 
data indicate that the results are applicable at any spanwise 
station. 

Camber effect .— To analyze the effect of the interfer- 
ence between the tunnel walls and the airfoil camber, the 
thickness and wake of the airfoil are considered to be removed 
and the airfoil reduced to its mean camber line. The result— 
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lug Infinitesimally thin airfoil may then he replaced by a 
sheet of continuously distributed, bound vortices which, In 
the general three-dimensional case, consist of both spanwise 
and chordwlse vortices. The velocity Induced at any given 
point on the camber line 1 b then obtained by Integration ovc-r 
the entire vortex sheet. Ae in all thin— airfoil theory the 
distribution of bound vorticity must be such that the result- 
ant of this induced velocity and the free— stream velocity is 
tangential to the camber line at all points. As will be seen, 
however, the actual theoretical determination of the distri- 
bution of vorticity is not necessary In this case. 

In calculating the velocity field of the vortex system, 

It is assumed that the bound vorticity is distributed in the 
middle plane of the tunnel — that is, in the xy— plane — 
rather than along the camber line and that the induced veloc- 
ity at any point on the camber line is the same as the induced 
velocity at the corresponding point in the xy— plane. Prom 
equation (5), the velocity potential at any point x, (u , 9 
far a vortex element on the y-axis at the point y = h (9 0 
= 0, u) 0 = tj ) is 



The term for m = 0 disappears by virtue of the factor m 
in the numerator of the general term. The vertical Induced 
velocity v z 5 in the incompressible fluid is then 

i _ dtp = dip d9 dtp dm 

8 ds dS dz duidz 


Por points in the xy— plane (0 = 0, <u 
Thus, at points in the xy— plane, 



0 and 

6z y 


T i - _ t 

Z a 

rrr s 


r* dn y V * a 6 B j ffl U B y) J a u a n) 

nr® i i /_i 


m=l 8=1 _ 'N * 2 T 3/i 

ny v - x7w 8 ” ( ' 


The complicated double series in this equation can be reduced 



HAG A ACE Ho, 5D31 


9 


to a single series and. the mathematics of the problem greatly 
s Amplified., "by . limit ing the discussion to the chordwise sec- 
tion of the airfoil at the center line of the tunnel (y = 0, 
z = 0). From the known relations for Bessel functions ( cf . 
reference 13 ) f for y = 0 

J i(*sy) _ . 


Jm(^-sy) „ 

= 0 for m > 1 


Thus, at points on the x— axiB , 


T i = _ r>liT 1 V e Xfl3C J^Xgn) 

z 2nr» L-. , X 

• 8=:L n x * J i a cv) 

8 

whers the summation with respect to s extends over all 
the positive roots of the equation 

J i *( ^a r ) “ 0 


From Bessel's differential equation 


J i"< x Br> = - (l- 


( 10 ) 


where the double prime denotes the second derivative of the 
Bessel function with respect to its argument. Equation (8) 
can thus be written 


r'dn y 9 XflX 

2nrq \ B r J 1 (A B r) j/'llgr) 


( 11 ) 


As mentioned, this equation is valid only for negative values 
of x. 
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It 1b apparent that tlie Berlea In equation (ll) is rapid- 
ly convergent for large negative values of x, but that the 
convergence is alow for small negative values. Since In the 
evaluation of the velocity induced by the vortex sheet it la 
the small values of x which are of primary importance, equar- 
tion (ll) cannot be applied directly lr. the present case. It 
is possible, however, by means of a method demonstrated by 
Vatson (reference 15), to express the series of this equation 
as a combination of elementary functions and a Berios of as- 
cending powers of x and ti . The resulting series is readi- 
ly applicable to the present problem. 


The detailed procedure for the transformation sf the 
series of equation (ll) 1 b given in appendix A. By applica- 
tion of the final result, equation (11; may be written 




A. B X 


*nrTi — 

t> 


_ r 1 dp 

2nrp 


ak+i 


rx 


2tvn/t 


p 3 +x £ 


=i X a T J i(*s r ) J i ^ X s r > 

* V y ( ~ 1}> 1 

k tr 0 kl ( k+l)'- (2p+l)i 2 ak+1 r sC k+ P+i) 


r z + p a 
2rp 

ap+i 


(12 


The double summation extends over all integral values of k 
and p from zero to positive infinity. The numerical coef- 
ficient ** 2(1 -+t>+ 1) = M> 2 f is e lven tha Integral 




2f 


1 

( 2f + 1 ) TT 



t sf 2 dt 


[i a '(t)3 


(13) 


Here I x (t) is a modified Bessel function of the first kind 

of order unity, and I x # ( t ) denotes the derivative of Ij(t) 

with respect to its argument. The numerical values of u.' . 

af 

for f = 1, 2, 3, 4 are evaluated by means of a series ex- 
pansion in appendix A. 


It is readily shown that the first term on the right- 
hand side of equation (12) agrees with the induced velocity 
computed for x = 0 by the more elementary theory of tunnel- 
wall interference which considers only the effects of the 
trailing vortices and their images. To this end, consider 
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the two-dimensional flow in a plane normal to the axlB of tho 
tunnel an Infinite distance downstream (fig. 2). The theory 
states that the 1 induced v elocity ■ at ■ a given point-(y,s) in 
this plane is twice as great as the Induced velocity at the 
corresponding point in the plane x s 0 (cf. reference 12, 
p. 260"). In the plane x = «» , the trailing vortices of the 
U— shape vortex previously considered constitute a vortex pair 
having an Infinitesimal spacing dT) and situated at the 
point y = T} , a » 0, The circulation of each vortex of the 
pair is T 1 and is directed as shown in figure 2. The bound- 
ary condition that there shall be no flow normal to the wall 
of the tunnel can be satisfied by the introduction at the 
point z = 0, y =* r a /Ti of an image vortex pair with a spac- 
ing r s dTi/Ti a and with the circulation of the vortices directed 
as Indicated. The vertical velocity induced at the center of 
the tunnel by the trailing vortex pair is 


Vst i ,(a>) 


r 1 dn 
2 ttti 3 


and the vertical velocity induced at the same point by the 
image vortex pair is 



The total vertloal velocity at the center of the tunnel at 
x - 00 is then the sum of these two velocities; that' is, 


r»dn (r 8 +ri a ) 
2tt r a n a 


(14) 


The vertical velocity at the center of the tunnel at x = 0 
is one— half of this value, or 


▼«'( 0 ) 


r»dT) (r a +T] a ) 

- 4 tt r a Ti a - 


(15) 


This value agrees with the result of equation (12) for the 
special case x = 0. Thus, the first term on the right-hand 
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side of equation (12) represents the vertical Induced velocity 
on the center line of the tunnel at x = 0 and 1 b attribut- 
able entirely to the trailing vortices and to the Interference 
between thess vortices and the tunnel walls. !The remaining 
terms represent the variation in Induced velocity due to a 
displacement a distance x upstream from the origin. These 
terms arise both from a change ih the effect of the trailing 
vortices and their wall interference and from the now— active 
effect of the transverse bound vortex and its Interference 
with the tunnel walls. 

Although equation (12) was deduced for negative values 
of x, it can be shown that it is applicable to positive 
values of x as well. According to Ton Karman and Burgers 
(reference 12, p. 2&7 ) , the vertical induced velocity at 
-x is related to the corresponding velocity at +x by the 
equation 


v z '(-x) = v z »(oo) - vj-'t+x) 

By virtue of this relation, together with the fact that 
v z '(0) = ^ t b 1 (°°) » it follows that 

v z '(+x) — v z '(C) = — [v a '(— x) — v z ’(0)] (16) 


That is, the difference between the induced velocity at a 
given station x and the induced velocity at x = 0 must 
be an odd function of x. The terms containing x in equa- 
tion (13), which were derived to represent this difference 
for negative values of the variable, .are seen to constitute 
precisely Buch a function. Thus the expansion of equation 
(12) is valid for positive as well as negative values of x. 

The vortex sheet which represents the entire airfoil can 
now be built up by the superposition of UV-shape vortices in 
the xy— plane, and the total Induced velocity found by inte- 
gration of equation (12) over the entire system. The lead- 
ing edge of the airfoil is placed on the y— axis as shown in 
figure 3; the trailing edge then lies at x = c, where c 
is the chord of the airfoil. The circulation of an elementa- 
ry vortex having an infinitesimal span dT) and situated at 
the point x = I, y = n is taken to be (dr’/d£)d|, where 
(dn/dt) is the vorticity per unit length of the chordwise 
section at the Btation y = tj . The vertical velocity induced 
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at the chordwise station x on the center line of the tunnel 
by a single elementary vortex la given hy equation (12) if x 
and t 1 '' are r eplaced hy (x— |) ’“and '( dT* /d£)d| ■, respectively. 
The total vertical velocity induced hy the complete airfoil 
is then given hy the double integral 



r®+n® r( x— £) 

TT| a T)®^® + ( x- £) a 


CO CO 


'll 

ir =n rv = f 


( - 1)P ■‘View) ^ ■ 

k! ( k+ 1 ) l ( 2p+l ) l 2 sk r a ( k+ P + i) 


dn d I 


(17) 


The integration of equation (17) requires a knowledge of 
(dT 1 /d() as a function of ti and (. Theoretically, 

(dT'/dO could he determined from the requirement that the 
induced vertical velocity at every point on the camber lino 
must bo auch that the resultant of this velocity and the 
free— stream velocity is tangential to the camber line. This 
method of procedure leads, however, to a complicated double 
integral equation, the splutlon of which does not appear 
feasible. Some assumption concerning the distribution of 
vorticity must therefore be made if the problem is to be 
solved. To aid in the cholco of n suitable assumption, ex- 
periments were carried out to determine the pressure distri- 
bution, both chordwise and. spanwlse, over an airfoil spanning 
a clcs ed— throat circular tunnel. The airfoil used in the 
experiments, which are described in detail later in the re- 
port, had an HACA 4412 section and was untwisted and of con- 
stant chord. The results of these experiments reveal that 
for such an arrangement the lift is sensibly uniform across 
the span for angles of attack below the stall. This fact 
is illustrated in figures 6 and 7, which show the experimen- 
tal spanwlse lift distribution for the airfoil at various 
angles of attack in wind tunnels affording chord— diameter 
ratios of 0.357 and 0.625. These results were at first re- 
garded as rather surprising. Later, however, it was realised 
that they are only what might logically be expected from gen- 
eral considerations of the conditions of flow in a closed- 
throat tunnel. A demonstration ©f this fact is given in ap- 
pendix B, in' which it is. shown that the. lift distribution is 
uniform across an untwisted, constant— chord airfoil spanning 
any closed— throat wind tunnel, irrespective of the cross-sec- 
tional shape of the tunnel. Detailed examination of the pres- 
sure distributions from which the results of figures 6 and 7 
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were obtained reveal further that at a given angle of attack 
the chordwise pressure distribution is sensibly the same for 
all spanvlse stations on the airfoil; that is, the lift per 
unit chord at any given chordwise station is constant across 
the span. It 1 b to be expocted that this result, though ob- 
tained for a particular airfoil, will be equally true for any 
ordinary camber— line shape. Thus it is reasonable to assume 
that the distribution of bound vortlcity is not a function of 
the spanvise position on the airfoil; that is (dT'/dl) is 
independent of t> . 


On the basis of this assumption, equation (17) may be 
wr itt en 


o — r 


r 2+ P a _ r( x— £ ) 

TTl 3 TlVn a +( Tr- £)S 


no co 


(~ p p ^ g(t+p+1) ( g - t )*** 1 

kit) p^6 ki (k+l)i ( 2p+l ) l 2 afc r a( k+ P +1 ) 


t y 


dtl d5 (18) 


and the integration carried out with respect to p . The first 
tws terms of the integrand, however, become infinite at the 
point n = 0. These singularities, which are due to the ef- 
fects of the vortices trailing from the vortex elements on 
the r-axis, require that special care be taken in the inte- 
gration. The evaluation of the integral must be carried out 
from — r to — e and from + e to +r , and to the resulting 
function must be added the effects of the trailing vortex 
pairs of span 3 c which straddle the x— axis . The limit of 
this sum must then be taken as e tends to zero. The verti- 
cal velocity Induced at the point x on the x— axis by the 
vortices trailing from a vortex element of span 2c symmetri- 
cally placed at x = £ , y = 0 Is 


v 


i 

z 


J- feVi + a.(*-g) 1 

4tt Vdg/Lc c v€® + ('x-|)"*"J 


d t 


Since the first two 
order powers of p , 
+c to +r will be 
with respect to p 


terms of equation (18) contain only second— 
the Integrals from — r to — e and from 
equal. The integral of these two terms 
thus becomes finally 
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The Integration with respect to T| of the double series in 
equation (18) presents no difficulty. The expression for 
v z ‘ thus becomes after integration 
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For constant spanvlse circulation, the trailing vortices 
finally disappear in the Integration with respect te T) . 

The integrand of equation (19) thus represents the Increment 
of vertical velocity Induced by an elementary vortex of con- 
stant circulation completely spanning the tunnel. 

It will now be assumed that the chord of the airfoil is 
small enough as compared with the dimensions of the wind 
tunnel that powers of ( x— 5)/r greater than the first may be 
neglected in the integrand of equation (19). This is equiva- 
lent to assuming that powers of the chord— diamet er ratio 
( c/d) higher than the second may be neglected in the final 
equations for the tunnol— wall corrections. The approximation 
is accomplished by expanding the first term of the integrand 
in ascending powers of (x— |)/r and discarding all terms con- 
taining powers higher than the first and by retaining only the 


ft > 
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= o terms of the double series. This gives for the induced 
elocity 


_ _i_ r dr* 

4ur J d{ 


(¥) 


1 + 


m 

w 

y ■ M» ^ g( k+ 1 ) 

^6 kl 


( k+ 1 ) ! ( 2k+l ) 2 


ak-i 


dt 


which may be written 



_1 + y >* ** 2( k + 1 ) I 

2 k z J 0 kl (k+l)» ( 2k+l ) 3 ak j 



By substituting the numerical values for the coefficients 
^ *2( k+1 ) from equations (A2L0) of appendix A, this equation 

may be written to an accuracy of three significant figures as 




m hf (so) 


The foregoing result, which was derived by assuming the 
fluid to be incompressible, can be modified for the effect 
of compressibility by the methods of Goldstein and Young. 

The modification is most readily performed by means of the 
so-called "Method II" (reference 4, pp . 5—6), which compares 
the compressible and Incompressible flows for equal values of 
circulation. If the Mach number of the compressible flow at 
the position of the airfoil is denoted by M, it is readily 
shown on the basis of this method that for a given distribu- 
tion of vorticity the vertical velocity induced in a compres- 
sible fluid at any p oint on the center line of a tunnel of 
radius r is J 1— M® times the corresponding velocity at the 
sam e p oint in an incompressible fluid in a tunnel of radius 
v/l— M a . Thus, from equation (20), the vertical velocity v z * 
in a compressible fluid in the actual tunnel of radius r is 
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-( 21 ) 


The first term of this equation represents the vertical veloc- 
ity that would he induced hy a vortex sheet of Infinite span 
in an unlimited fluid field. Ihe second term thus represents 
the Interference effect of the tunnel wall. 


Equation (21) may 
suit from reference 5, 
for an airfoil in a cl 
tunnel. After alterat 
sign conventions of th 
ence 5 gives for the v 
an infinitesimally thi 
a t wo— dimens ional— f low 


he compared with the corresponding re— 
which discusses the wall interference 
osed— throat two— dimensional— flow wind 
ion to conform with the notation and 
e present paper, equation (41) of refer- 
ertlcal velocity at the camher line of 
n airfoil mounted on the center line of 
tunnel of height h 


= 


-m 3 f /ar»\ r i _ 

! rr -J \d £ ) L X— t 


6h a ( 1— M a ) 




( 22 ) 


Comparison of equations (2l) and (22) shows that an infinites- 
imally thin airfoil spanning a clos ed— throat circular tunnel 
of radius r experiences at its mldspan section the same in- 
terference as would he experienced hy the same airfoil in a 
cl os ed— thr oat two— dimens ional— flow tunnel of height 


hn - 


,*/6( 0.579 


r = 1.686 r 


or, in termB of the tunnel diameter d, 


h x = 0.843 


(23) 


This result makes the later determination of the interference 
corrections for the circular tunnel very simple, since the 
corrections for the rectangular tunnel are already known. 

It is readily shown hy means of equation (6) that the 
vortex system which represents the infinitesimally thin air— 
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foil induces no axial velocity at any point in the xy— plane. 
It follows that airfoil camber has no effect upon the axial 
velocity or pressure gradient at the position of the model. 


Thickness effect .— The interference effects associated 

with airfoil thickness can be found by reducing the given 
airfoil to its base profile and analyzing the interference 
between the tunnel wall and this profile. The base profile 
is defined as the profile the airfoil would have if the cam- 
ber were removed and the resulting airfoil placed at zero 
angle of attack. If it is assumed that no wake is present, 
the interference between the tunnel wall and this symmetri- 
cal airfoil can be found by applying the results ef equation 
(5) to Lock's method of analysis of the interference on a 
symmetrical body in two-dimensional incompressible flow. 
(Lock's original analysis appears in reference 1; an alterna- 
tive explanation of the method is given by Glauert in refer- 
ence 2, pp . 52—57.) 

• 

Lock's method of analysis, which assumes that the chord 
of the airfoil 1 b small as compared with the dimensions of 
the tunnel, consists essentially in replacing the given 
symmetrical airfoil by an equivalent two-dimensional Bource— 
sink dotiblet and calculating the interference between thiB 
dbublot and the tunnel boundaries. The strength of the doub- 
let in any given case is proportioned so that it induces at 
a considerable distance from itself in free air a velocity 
equal to the velocity induced at the same point by the origi- 
nal airfoil. In the two-dimensional case, the Interference 
flow at the position of the airfoil is then readily found by 
introducing an infinite series of images of the doublet such 
as to satisfy the condition that there shall be no flow nor- 
mal to the tunnel boundaries and calculating the velocity 
induced at the airfoil by this system of images. For an air- 
foil spanning a clos ed— throat rectangular tunnel at mid- 
height, the net result of the wall Interference for the in- 
compressible case is to Increase the effective axial velocity 
at the position of the airfoil by the amount 




TTH 

6h a 


where ^ is the strength of the doublet used to represent 
the airfoil. It is shown in references 4 and 5 that the ef- 
fect of fluid compressibility is to increase this interference 
velocity by the factor l/[l— (M') a ] s/s, where M' is the 
Mach number of the undisturbed stream in the tunnel. Thus, in 
the compressible case, 
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jni 

6h a [l-(M«) S ] 5 / 8 


( 24) 


The problem of the symmetrical airfoil in a closed- 
throat circular tunnel can also he solved by replacing the 
airfoil by an equivalent doublet spanning the tunnel. In 
this case, though, the Interference for the doublet cannot 
be found by the method of images. If the doublet UBed is 
composed, however, of two vortices In a plane normal to the 
stream Instead of the customary source and sink In line with 
the stream, the interference velocity oan be calculated by 
means of equation (5). Since the velocity fields of the two 
types of doublets are identical, the interference calculated 
by means of the vortex doublet Is the same as that which 
would be obtained if the source— sink doublet were used. 


Consider a vortex element of circulation T* and Bpan 
dr) at the point u> 0 , 0 O in the yz— plane (fig. l). From 
equation (5), the streamwise velocity v z ’ induced at any 
point x, to, 6 upstream from the origin by this element and 
its accompanying trailing vortices is 


Vj 1 = 


dtp 

dx 


r»dri 

nr a 


e J n (k s u)) 


00 . 00 

y y - 

m=0 ^1 (l- — Hi 


r j_( ) 

m cos 0 O sin m( 0— 0 0 ) 

L u > 0 

+ \ a sin 0 0 cos m(0-0 Q ) J n '(^ B (o 0 )J 


(26) 


At a point in the middle plane of the tunnel ( xy—plane) , 
0 *= 0, to *= y, and the velocity is 


v x' - 


£ 

nr 


GO CO 


e Jm( ^ s7 ) 


m=08 ^ 1 C 1 - 

X.j^- m cos 0 o sin m0 o — + A 8 a in 0. 0 cos m0 o Jm^C^s^o)^] 


( 26 ) 
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As before, the double series In this equation reduces 
tc a single series if the discussion is limited to the inter- 
ference at the center line of the tunnel. For points on the 
center line, 


JoU 8 y> = J o(°) 55 1 

JmUey) = J m (°) = 0 for m ^ 1 


and the streamwise induced Telocity becomes 

00 X x 

r'd« V 1 p 

v. 


I 


= _ r'dn 

3rrr r 


I 

8 = 1 


sin e 0 J o '(\ a u) 0 ) 


J 0 a (X 8 r) 


(27) . 


From the known relations for the Bessel functions 


J o’Us u, o) = ~ J i(*«“o) 


so that equation (27) may be written 

„ i = nm y 


e X ® x sin 6 o J l ( x 8 (u o^ 


2nr ' 


8 = 1 


Jo S ( x «r) 


(28) 


(29) 


As required by equation (4), the summation with respect to 
s in this equation extends over all positive roots of the 
equat ion 


Jo'( X B r > = “ J l( X 8 r) ■ 0 


(30) 


As the next step, consider a pair of symmetrically placed 
elementary vortices composed of a vortex of circulation — r* 
at the point u> 0 , 8 0 and a vortex of circulation +r 1 at 

the point u> 0 , — tf 0 . From equation (29), the streamwise veloc- 

ity induced at a point on the center line of the tunnel by 
this vortex pair and the accompanying trailing vortices is 

, _ _ r'dn V 1 eXg3C sln e o J i( X B<»o) 

Z TTT 3 L, J 0 2 ( X S r ) 

8 = 1 


( 31 ) 
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which may also he written 


2H u) o ain 9. 0 dt) y e 8 J^^a^o) ( 32 j 

Sur® °>o 8 4 i Jo a U a r) 

The expression (2r ( u> 0 Bin 0 O ) which appears in thiB equa- 
tion is the product of the vortex strength and the distance 
between the vortices. 

How let the distance between the vortices tend to zero 
while the vortex strength increases in such a way that the 
product (ET'ujQsin 0 O ) retains a constant value 4 . The re- 
sult in the limit is an elementary vortex doublet of strength 
and span dr) at the point ro 0 =* ti on the y— axis . The 
streamwise velocity induced on the-oenter line of the tunnel 
by this elementary spanwise doublet and the accompanying 
trailing vortex doublets is then 



m. y e XflI J^sh) (33) 

2 rrr s n ^ J o B ( X s r ) 

8 = 1 

As before, the infinite series in this equation is 
rapidly convergent for large negative values of x, but the 
convergence is slow for small negative values and is non- 
existent when x = 0. Once again, however, the series can 
be expressed as a combination of elementary functions and a 
power series which is readily applied to the problem at hand. 
The details of the transformation are given in appendix A. 

By means of the final result, equation (33) can be written 



Vj" = - 


1* 


2 nr a 
(-1) P 


3U a +x s ) 3 / a 


** b( k+p+i ) 


n ak x»P 


2 1 - 

k top=o ‘ti.OwDi (. 2p)i a 81 "-. 1 r ak+ >P +1 J 


(34) 


where the double summation extends over all Integral values 
of k and p from zero to infinity. The coefficient 
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^ a ( k+ p+l) = W a f is given "by the integral 

U , = 1 f ** 

af ( 3f+l ) n J I^tt) 

o 

The numerical values of this integral for f 
are evaluated in appendix A. 


(35) 

1. 2, 3, 4 


The induced velocity for a doublet spanning the tunnel 
is now readily found by taking the doublet strength p con- 
stant across the span and integrating equation (34) with 
respect to tj from — r • to +r . This gives finally 


v 


i 

x 


M; 

Srtr 2 


x 2 + x 2 


v~i ( l)^* ^a(k+p+i) xS ^ 

1~ 0 “ 0 k '* ( k +l)i (2p)». (3k+l) 2 sk r a P 


(36) 


In the integration across the tunnel, all the trailing vorti- 
ces, of course, disappear. 

It is apparent from the symmetry of the problem, that 
the streamwise velocity Induced by a doublet spanning the 
tunnel must be an even function of the variable x. Equation 
(36), which was derived for negative values of x, is seen 
to be such a function and is thus applicable to positive val- 
ues of the variable as well. 


The values of v x ' for vanishingly small values of x, 
that is, at the position of the doublet, is then found from 
equation (36) by expanding the first term in ascending powers 
of x/r and discarding all terms containing second powers 
and higher and by retaining only the p = 0 terms of the 
double series. This gives 


▼x 


i 



1_ _ y ^a( k+i) 

2 k Tb k * ( k+ l)! (2k+l) 2 sk 


(37) 


After substitution of the pumerical values for the coefficients 
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^s(k+i) from equations (ASS) -of appendix A, this equation 
'becomes to an accuracy of three significant figures 


r . = LL_ + 1-356 y, 

X 2ttx® 2rrr a 


(38) 


The first term of equation (38) is the velocity induced 
by a doublet of infinite span in a field of unlimited extent. 
The remaining term therefore represents the effect of inter- 
ference between the doublet and the tunnel vail. Thus the 
net result of the interference between the airfoil thickness 
and the tunnel wall for the incompressible fluid is to in- 
crease the effective stream velocity at the position of the 
airfoil by the amount 


a y. = lij 5 6. ja 

1 2nr a 


(39) 


In any particular case, p. is again equal te the strength 
of the doublet used to represent the given airfoil. 

The result of equation (39) can be modified for the ef- 
fect of fluid compressibility by the method *f reference 4. 

In this case, the modification is most conveniently performed 
by means of Method I (reference 4, pp. 3—5) which compares 
the compressible and incompressible flows for a given airfoil 
of unaltered shape and else. By this method, it is readily 
shown that the streamwise velocity Induced in the incompressi- 
ble fluid at any point on the center line of a tunnel of radius 
r is l/^/l“( M 1 ) a times the corresponding velocity at the 
same point in an incompressible fluid in a tunnel of radius 

iV 1— ( K ’ ) 3 • Here M’ is, as before, the Mach number in the 
undisturbed stream. The increment in- axial velocity in the 
compressible case is thuB 


A 1 7' 


1.356 u 

2nr a [l-(H')^ 3 / a 


(40) 


. . Comparison of equations (24) and (40) shows that, if no 

wake is present, a symmetrical airfoil spanning a closed- 
throat circular tunnel of radius r experiences at its mid- 
span section the same increase in axial velocity as would be 
experienced by the same airfoil in a clos ed— throat two— dimen- 
sional— flow tunnel of height 
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h a «= 0.779 4 (41) 

The foregoing result greatly simplifies the determination of 
the true stream conditions at the position of the airfoil in 
the circular tunnel since the necessary equations for the 
rectangular tunnel are already known. 

Consideration of the symmetry of the system formed by 
a base profile spanning the middle of a circular tunnel indi- 
cates that the Interference between the wall and the airfoil 
thickness does not Influence the vertical induced velocity 
v B ' at any point on the airfoil. Similarly, the airfoil 
thickness has no effect upon the streamwise pressure gradient 
in the tunnel at the position of the airfoil. 

Wake effect .— It is shown in general terms in reference 

5 that the interference between the wake of a body and the 
walls of a closed— throat wind tunnel gives rise at the posi- 
tion of the body to a velocity increment and a streamwise 
pressure gradient which are not present in free air. This 
is true for any type of body and any shape of tunnel test 
section. The magnitude of this velocity increment and pres- 
sure gradient in the case of an airfoil spanning a closed- 
throat rectangular tunnel can be determined approximately by 
replacing the wake by the flow from a suitable fluid source 
and the tunnel walls by an infinite system of image sources. 
In the case of the airfoil spanning a clos ed— thr oat circular 
tunnel, this treatment is no longer possible since no system 
of images sources is known which will satisfy the boundary 
conditions at the tunnel wall. A more complex method of 
analysis could conceivably be devised for this case; however, 
since the calculation is highly approximate even in the case 
of two-dimensional flow, such a.i analysis does not appear 
warranted. Tor present purposes it is probably sufficient to 
assume that the midspan section of the airfoil in the circu- 
lar tunnel experiences the same velocity increment and pres- 
sure gradient as a result of the wake Interference as does 
the same airfoil in a rectangular tunnel of a height defined 
by equation (4l). This assumption leads to the simplest ex- 
pression for the final correction to the measured drag coef— 
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ficlent and should gi re results yhich are reasonably accurate. 
If It is assumed that tha oenter of the wake lies in a hori- 
zontal plane containing - the diameter of the tunnel, it follows 
from considerations of symmetry that the wake interference 
does not contribute to the vertical induced velocity v z ' at 
the airfoil. 


It has already been indicated that the Interference asso- 
ciated with the camber cf the airfoil has no effect upon the 
stream velocity at the model. The total Increase in velocity 
for the complete airfoil in the oircular tunnel is thus given 
by the sum of the increments causod by the thickness and the 
wake of the airfoil. In reference 5 it is shown that for the 
analagous case of the airfoil in the rectangular tunnel, the 
true velocity V at the position of the airfoil may finally 
be written 


7 


7 « 



± Acr + 

C 1 — ( M • ) a ] 3 / a 


1+0.4(M» ) a 
l-( M » ) a 



(43) 


where a and T are factors dependent upon the size of the 
airfoil relative to the tunnel, A is a factor dependent 
upon the shape of the base profile, and C& 1 is the drag 

coefficient of the airfoil as measured in the tunnel. The 
first correction term in this equation represents the veloc- 
ity Increment caused by the airfoil thickness and is found 
by substituting the proper value for the equivalent doublet 
strength in equation (24). The second correction term repre- 
sents the velocity increment associated with the wake of the 
airfoil. 


The factors a and T in equation (43) are defined by 


and 


n a ( c 'N* 

(43) 

CT = M \hJ 

t ■*(■£■) 

(44) 


where (c/h) is the ratio of the • airfoil .chord to the tunnel 
height. An analytic expression for A Is given in equation 
(3) of reference 5. Values of A for a number of base pro- 
files are given in table I, which is reproduced from this 
reference. 
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If it ia aBBuned that the height of the equivalent rec- 
tangular tunnel with regal'd to the wake interference ie the 
same as that given by equation (4l) for the thickneBB inter- 
ference! the true velocity in the circular tunnel ie found 
simply "by subs t itut ing h a from equation (4l) for h in 


the factors o and T of equation (42). The 
at the midspan section of an airfoil spanning 


tunnel 

is 

thus 

> 

ii 

> 

{i 

+ 


L 

Cl— ( M 

wher e 

the 

factors 


) 3 

T * 


true velocity 
a circular 



(45) 


and 



( 46‘) 


(47) 


A correction to the stream velocity Implies corrections 
also to the stream dynamic pressure, Reynolds number, and 
Mach number. These corrections for an airfoil spanning a 
rectangular tunnel ha^e been determined in reference 5 on the 
basis of the assumption that the flow is adiabatic. The 
corresponding corrections for the circular tunnel can be found 
by replacing the factors T and a in equations (29), (32), 
and (33) of reference 5 by the factors T a and ct 8 of the 

present paper. The true dynamic pressure q, Reynolds num- 
ber R, and Mach number M at the midspan section in the 
circular tunnel are thus related to the correspond ing quanti- 
ties in the undisturbed stream (denoted by primes) by the 
equat i ons 


— Mm,) * | U8) 

L Cl-(M‘) a ] /* 1-(M*) S 2 J 

= » ( 1+ Ar-° •*■">* Aa. ♦ , 0d .\ ( 49 ) 

1 Cl-tM 1 ) a 3 3 / a 5 1 — ( H 1 ) * ’ f 

■ «' {i+ T ~£rn *>.+ ■ [1,MlM ' )i][lt0 ; 4( ”' ^ u°d'} (so) 

L Cl— ( M 1 J 5 ] 3 / a l - (M> ) a J 
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Burner ical values ctf the functions of M 4 which appear In 
these equations are given in table II, which Is roproduced 
from reference 5. " 

At low Mach numbers! the terms containing T a ca' in the 

equations for the corrected stream characteristics are usually 
negligible as compared with the terms containing A a a . At 

high Mach numbers, however, where the drag coefficient is very 
large, the terms with T^a* predominate. 


Relations between Airfoil Characteristics 
in Tunnel and in Pree Air 

The characteristics of the airfoil in free air are now 
readily determined in terms of the characteristics at the 
midBpan section in the tunnel. It is simply necessary to ap- 
ply the results of the preceding sections to the relatione 
already derived in reference 5 for the airfoil spanning a 
rectangular tunnel. 

Briefly, the method of reference 5 relates the section 
characteristics in the tunnel at an undisturbed stream veloc- 
ity V 1 to the characteristics in an unconfined stream hav- 
ing a velocity equal to the true velocity V which exists 
at the position of the airfoil in the tunnel. The relation 
is obtained on the basis of equal valueB of the so-called 
cotangent cemponent of lift in the tunnel and in free air, 
this being necessary to assuro that the essential character 
of the pressure distribution over the airfoil is the same in 
both cases. By this procedure corrections are derived which 
may be applied to simultaneously measured lift, moment, and 
drag coefficients and angle of attack in the tunnel to obtain 
the corresponding quantities in free air. These corrections 
appear as functions of the factors A and o, of the prod- 
uct Tea* , and of the Mach number M 1 of the undisturbed 
stream. The correction to the angle of attack, which arises 
out of the Interference effects associated with camber, is 
proportional to a and independent of A and Tea’. The 

correction equations for the lift and moment coefficients 
contain corresponding terms proportional to ct alone, to- 
gether with terms which depend upon the thickness and the wake 
effects and are proportional to the products Act and Tea ■ 
The correction to the drag coefficient appears as two terms. 


38 


HACA ACH Ho. 5D31 


proportional to A o and Tea' , respectively. The term pro- 
portional tn Aa is in this case composed basically of two 
parts, sne due to the thickness effect and one due to the 
wake effect. 


The correction equations for the airfoil spanning a 
circular tunnel can he derived directly by modifying the equa- 
tions of reference 5 in accordance with the results of the 
preceding sections. Since the terms containing a exclusive 
of A appear as a result sf the camber effect, the tunnel 
height h in such terms must be replaced by 0.843d as re- 
quired by equation (33). In the terms which depend upon the 
thickness and wake effects and are distinguished by the prod- 
ucts A a and tc^ 1 , the quantity h is replaced by 0.779d 
in accordance with equation (4l). This involves the assump- 
tion already mentioned that the height of the equivalent rec- 
tangular tunnel with regard to the wake effect is the same 
as that calculated for the thickness effect. 


As in reference 5, the free— air lift, quarter— chord- 
moment, and drag coefficients referred to the true dynamic 
pressure q are denoted by the conventional symbols. The 
corresponding quantities measured in the tunnel and referred 
to the apparent dynamic pressure q’ are denoted by the same 
symbols with primes added. The final equations for the cor- 
rected aerodynamic coefficients are then 
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(51) 


[ 3— ( M 1 ) g ] [1+0. 4( M 1 ) a ] T ,1 
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°d = 



3 -0.6(K' )* 

[l-(M« )*] 3 /» 


_ [2-( H 1 ) 8 ][ l+0.4( M * ) 3 ] 
1 - (M') a 



(53) 


and the corr eeponding angle of- attack In degrees is 


, 57.3(7, f , \ 

a = a* + ========= "S o j + 4 Cn / 


where the factor cr 1 le given by 

= 0-289 (jJ 


(54) 


(55) 


and the factors t 2 and a g are as already defined In 

eruptions (46) and (47). Numerical values of the compressi- 
bility factors which appear In these equations are given In 
table II. Tho corrected quantities correspond to the true 
Heynolds number and true Kach number as given by equations 
( 49 ) and ( 50 ) . 


From a rigorous standpoint, the foregoing corrections 
apply only to data obtained from chordwlse pressure distribu- 
tions at the uldspan section of the airfoil. Actually, as 
has already been pointed out in the discussion of camber ef- 
fect, the experimental chordwlse pressure distribution at any 
given angle of attack is sensibly oonstant across the span. 

The corrections should therefore be applicable with suffi— - 
'cient accuracy to data obtained from pressure distributions 
at any Bpanwlse station. 

Reference 5 also includes a method for correcting experi- 
mental chordwlse pressure distributions to free— air conditions 
in the case of an airfoil spanning a rectangular tunnel. The 
same procedure may be applied to -pressure, dis.tr lbut ions over 
an airfoil spanning a circular tunnel if the factor t is 
replaced by T a and the factor a by ct x wherever it ap- 
pears alone and by or s where it appears in the product Acr. 
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Choking at High Speeds 

As explained in reference 5, for tests of a model in 
any clos ed— thr oat wind tunnel, there is some value of the 
Mach number M' of the undisturbed stream which cannot be 
exceeded irrespective of the power input to the tunnel. This 
follows from the fact that at high speeds the combination of 
model and wind tunnel acts essentially as a converging— di- 
verging nozzle, and the flow in the tunnel exhibits the char- 
acteristics of the flow in such a nozzle. Thus, at some Mach 
number less than unity in the undisturbed stream, sonic veloc- 
ity is attained at all points across a section of the tunnel, 
usually in the vicinity of the model. When this occurs, in- 
creased power input to the tunnel serves merely to extend the 
region of supersonic flow dewnstream of this sonic section 
and haB no effect upon the velocity of the stream ahead of the 
model. The tunnel is then said to be "choked," and the Mach 
number M' of the undisturbed flow ahead of the airfoil has 
its maximum attainable value. This value is described as the 
apparent choking Mach number, the word "apparent" being used 
to differentiate this value from the corresponding free— air 
Mach number M which would be computed from equation (50). 

If it iB assumed that the section of sonic velocity is 
coincident with the section of minimum area between the model 
and the tunnel walls, the apparent choking Mach number can be 
obtained from olementary considerations of unidimensional 
adiabatic flow, as shown in reference 5. Tor the present case 
of a cons tant— chord airfoil spanning a circular tunnel, the 
apparent choking Mach number M’ ^ if finally defined for 

air ( *v = 1.4) by the relation 



where te is the "effective" thickness of the airfoil and 
d is , as before, the diameter of the tunnel. A graph of this 
relation is given in figure 4. As a matter of interest, the 
results are shown for the supersonic as well as the subsonic 
flow regime. The region above the curve represents an im- 
possible state of flow. 

In estimating the apparent choking Mach number in any 
practical case it is necessary to replace the effective thick— 
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neas t e by the projected, thickhess tp of the airfoil nor- 
mal to the direction of flow. Ae indicated in reference 5. 
this procedure leads, in the case of the subsonic wind tunnel, 
to an over es t i'mat ion of because it neglects the possi- 

ble contraction ef a portion of the Btream aft of the airfoil 
as well as the effect of the airfoil boundary layer. 

She importance of the boundary layer and the accompany- 
ing drag with regard to tunnel choking is pointed out in ref- 
erence 5, where the apparent choking Mach number is calculated 
for a flat plate at zero angle of attack in a two-dimensional- 
flow wind tunnel. Since the projected thickness for the plate 
is zero, the unidimensional theory would Indicate that no 
choking occurs. Actually, beoause of the faot that the plate 
experiences drag, choking does take place. Similar considera- 
tions hold, of course, for a flat plate spanning a circular 
tunnel. In this case the apparent choking Mach number for 
air (Y = 1.4) is given by the equation 


i . (<l 

tt \a 



) 


drag on choking for supersonic as well as subsonic wind tunnels 
is shown. It can be demonstrated that the points on the curve 
correspond to a Mach number of unity in the flow far down- 
stream of the model where the wake has spread completely to 
the tunnel wall. PointB above the curve represent impossible 
conditions of flow. In most cases encountered in subsonic 
tunnels, the apparent choking Mach number determined by the 
thickness of the airfoil and defined by equation (56) is usu- 
ally the lower. For very thin airfoils at small angles of 
attack, however, the value of M 1 c h given by equation (57) 


can have the lower value. At present no' way is known to com- 
bine the thickness and drag effects in a single calculation 
as should logically be done. > 


It should be noted, as pointed out in detail in refer- 
ence 5, that the flow in a tunnel at choking does not corre- 
spond to any flow in free air. Furthermore, for a range of 
Mach numbers just below choking, where the flow is influenced 
to any extent by the restrictions which finally promote chok- 
ing, any wall— interference correction is of doubtful accuracy. 
This is particularly true if the model is at an appreciable 
angle of attack so that sonic velocity is attained across the 
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atream on one side of the airfoil before It le on the other. 


EXPERIMENT 


The experimental invest igat ion was initiated for two 
reasons: (l) to determine the spanwise distribution of lift 

over an airfoil spanning a clos ed— thr oat circular tunnel, and 
(2) to examine the validity of the theoretical interference 
corrections derived in the preceding analysis. As has been 
previously mentioned, the development of the theoretical re- 
lations requires a knowledge of the variation in lift over 
the span of the airfoil. Since no theoretical or experimental 
evidence regarding this matter was available, the spanwise 
variation in lift was investigated experimentally for an NACA 
4412 airfoil for two ratios of airfoil chord to tunnel diam- 
eter. The results of these tests are also directly applicable 
to the examination of the validity of the theoretical cor- 
rection equations. 

The experimental work was performed in a low— turbulence, 
nonreturn— type wind tunnel with interchangeable throat sec- 
tions of 14— and 8— inch diameter. The two chord— diameter 
ratios were obtained by testing the same airfoil in each 
throat section. Since the airspeed was held constant through- 
out the tests, this arrangement permitted the Reynolds number 
and the Mach number to be duplicated simultaneously for the 
two chord— diamet er ratios. In this manner the effects of 
any variation in these parameters were eliminated from the 
tests. 

The NACA 4412 airfoil was used because a model of suit- 
able size was already available Ideally equipped for pressure- 
distribution tests. The model, which is described in refer- 
ence 16, was of 5— inch chord and 30— Inch span. This cherd, 
together with the two throat diameters, gave chord— diameter 
ratios of 0.357 and 0.625. In the tests, the airfoil ex- 
tended through the walls of the tunnel and was clamped in 
tight-fitting support blocks which prevented any leakage of 
air at the walls. The 54 pressure orifices located around 
the surface of the midspan section of the model were connected 
to a mult ipl e— tube manometer for measurement of the pressure 
distribution over the airfoil. To secure as accurate pres- 
sure-distribution data as possible, alcohol was used as the 
manometer fluid and the liquid heights were recorded photo- 
graphically. 
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Freseur a— distrlbut ion records vere secured at each of 
eight angles of attack from —4° to +15° at a Reynolds number 
of. ■ approximately 450,000 and a Mach number of approximately 
0.2 with the model mounted in bath the 14— inch and the 8— 
inch diameter throats. ■ The spanwlse distribution of lift 
was determined for each angle of attack by sliding the pres- 
sure orifices laterally from one wall to the other and re- 
cording the Indicated pressure distributions at a number ef 
spanwlse stations. The chordwlse pressure distributions were 
plotted and mechanically integrated to obtain lift and quar- 
ter-chord moment coefficients. Ho drag coefficients were ob- 
tained because the experimental installation did not permit 
balance measurements to bo made and wake surveys were not 
f eas lble . 

By testing the airfoil in both erect and inverted atti- 
tudes the inclination of the air stream with respect to the 
tunnel axis was determined for each throat section. The 
stream angle was found to be +0.45° for the 14— inch throat 
and 0° for the 8— inch throat. Corrections have been applied 
to all angles of attack for the measured angularity. 

The spanwlse distribution of lift coefficient uncorrected 
f 6r tunnel— wall interference is shown for the two chord— diam- 
eter ratios in figures 5 and 7 in which lift coefficients at 
various angles of attack are plotted as a function of the 
spanwlse location of the measurement plane. 

Curves of lift coefficient against angle of attack for 
the two chord— diameter ratios are shown uncorrected for tunnel- 
wall Interference in figure 8(a). The results given pertain 
to the section of the airfoil at the center line of the tunnel. 
The corresponding curves corrected for wall interference by 
means of equations (51) and (54) are shown in figure 8(b). 

In applying the ■ correct ions , the term containing Tea.' was 
necessarily omitted as no measurements of drag were made. 

For the values of c^ 1 to be expected in such tests, however, 
this term would be negligible in comparison with the remain- 
ing terms so that this omission does not affect the final re- 
sults. For purposes of comparison, section lift characteris- 
tics as obtained by Pinkerton from tests of a finite— span 
rectangular airfoil in the Langley var iable— dens ity wind tunnel 
(reference 17) are. .also shown.. These data correspond to an 
effective Reynolds number of 450,000 and are thus directly 
comparable to the results of the present test. 

In figure 9(a) curves of quar t er— chord moment coefficient 
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against lift coefficient are shown unesrrected for tunnel— 
wall Interference fer both- chord— diameter ratios. The same 
data are plotted in figure 9(h) after correction for wall in- 
terference by means of equations (5l) and (52). Also shown 
for comparison are the corresponding data from reference 1*7. 


DISCUSSION 


An examination of figures 6 and 7 reveals the previously 
mentioned fact that there is n.o appreciable variation in lift 
over the span of the alrfeil at all angles of attack up to 
those closely approaching the stalling angle. This observa- 
tion holds for both chord— diamet er ratios. In the vicinity 
of the stall a spanwise variation in lift appears which be- 
comes progressively mere erratic as the angle of attack is 
Increased. As might be expected, this variation becomes ap- 
parent at a lower angle in the case of the larger chord— diam- 
eter ratio. The results of figures 6 and 7 corroborate the 
conclusion of appendix B for the particular case of the air- 
foil spanning a circular tunnel. 

Prom figure 8(b), it is seen that the corrected lift 
curves for the two chord— d iamet er ratios agree almost exactly 
with one another except at angles near the stall. Below the 
vicinity of the stall the corrected data coincide with the 
results of reference 17 except for a constant angular dis- 
placement of approximately 0.2°. In reference 16, Pinkerton 
estimates that his values for the angles of attack may be 
too large by a constant error of appr sximat ely 0.25° because 
of a possible error in the assumed direction of the stream. 

It is thought that the angles of attack of the present ex- 
periments are accurate to within ±0.1°. These limits of ac- 
curacy are sufficient to account completely for the apparent 
angular displacement. 

In the region of the stall, the corrected lift curves 
for the two chord— diameter ratios do not mutually coincide, 
but the data for the chord— d iamet er ratio of 0.357 agree with 
Pinkerton's results within 2 percent. As previously men- 
tioned, Pinkerton's tests were made with a finite— span rec- 
tangular airfoil, for which the cross— span variation in lift 
is necesaar ily large. It is not to be expected that the de- 
termination of maximum section lift from such tests would be 
as accurate as from tests nf a through model, for which the 
cross— span lift variation is small. 
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It 1b seen from figure 9(b) that the corrected moment 
curves agree satisfactorily with each other and with the re— 
suits of r af er ence - 17 . 

In summary, for angles of attack helov those in the 
region of maximum lift, the results presented in figures 8 
and 9 demonstrate the validity of the theoretical lift, mo- 
ment, and angle— of— attack corrections for low Mach numbers 
and chord— diamet er ratios up to at least 0.625. Jor angles 
In the vicinity of maximum lift, the corrections are not 
strictly applicable up to such a large chord— diam et ar ratio. 
The results of the present test indicate that an accurate de- 
termination of maximum lift can be made with a chord— diameter 
ratio at least ae high as 0.36. An evaluation of the accu- 
racy of the correction equations at high Mach numbers is not 
possible on the basis of the experimental evidence available 
at present. It is to be expected, however, that the maximum 
permissible chord— diameter ratios will decrease as the Mach 
number incr-e&ses. 

The data of the present paper enable no definite con- 
clusions to be drawn regarding the validity of the drag cor- 
rection. However, in view of the accuracy of the other cor- 
rections for the circular tunnel and in view of the fact that 
the corresponding drag correction for a two— dimens ienal tunnel 
is known to be accurate, it is to be expected that this cor- 
rection will give a satisfactory evaluation of the wall inter- 
ference upon the measured drag. 

The equations of the present paper should not be expected 
to give accurate results when applied to tests in whioh air 
leakage occurs at the tunnel walls. In such tests the lift 
at the wallB drops markedly, bo that the assumption that the 
lift is uniform across the span is no longer valid. The im- 
portance of avoiding such leakage, if reliable airfoil char- 
acteristics are to be obtained, is pointed out in reference 
5 with regard to tests in two-dimensional tunnels. The same 
general considerations apply In the case of an airfoil span- 
ning a circular tunnel. 


CONCLUSIONS 


Airfoil data obtained from teste at subsonic speeds of 
an airfoil spanning the center of a closed— throat circular 
wind tunnel can be corrected to free— air conditions by means 
of the following equations; 
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where t s , a * , and ct_ are given bj 
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“ (f) 

(46) 

- °- 289 (f) a 

(56) 

«= 0.339 

(47) 


and A le & dimens lonlesB factor the value of -which depends vcpcn 
the ehape of the base profile of the airfoil. (See table I 
and equation (3) of reference 5.) The remaining symbols are 
defined in appendix C. Numerical values of the functions of 
M 1 which appear in these equations are given in table II. 
Experimental pressure die tr lbut Iona can also be corrected by 
proper modification of the method of reference 5 as Indicated 
in the text. 


Tests of an NAGA 4412 airfoil at low speed for two ratios 

of airfoil chord to tunnel diameter demonstrate the validity 

of the foregoing equations at low Mach numbers. At angles of 

attack below the region of maximum lift, the equations are 

applicable for chord diameter ratios up to at least 0.625, the 
maximum ratio tested. In the region of maximum lift a chord— 
diameter ratio of 0.35 is known to be permissible, and still 
higher ratios may give satisfactory results. An examination 
of the validity nf the equations at high Mach numbers is not 
possible at present, but the maximum permissible chord— diam- 
eter ratios may be expected to decrease as the Mach number 
lncr eas es . 

The tests also indicate that at low Mach numbers the 
spanwise lift distribution on an airfoil spanning a closed- 
throat circular tunnel is essentially constant except at 
angles of attack in the immediate vicinity of the stall. 

This result corroborates the general conclusion of appendix 
B, in which it is demonstrated that the lift is uniform across 
an untwisted, constant— chord airfoil spanning any closed- 
throat wind tunnel, irrespective of the cross-sectional shape 
of the tunnel. 

The correction equations cannot be expected to apply at 
or in the immediate vicinity of the choking Mach number, which 
is the maximum Mach number -attainable with a given combina- 
tion of airfoil and tunnel test section. The choking Mach 
number can be estimated by means of equations given in the 
report . 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field , Calif, 
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APPENDIX A 


TEASSPCBMATION OP SEBIIiS OP BSSSBL FUNCTION 


The series involving Bessel functions which appear In 
the discussions of the interference effects associated with 
airfoil camber and thiokness are, as pointed out in the text, 
poorly suited for use at small values of the variable x. It 
will be shown here, by means of a method demonstrated by 
Watson (reference 15), that the series may each be exuressed 
as a combination of elementary functions and a convergent 
power series. The resulting series are well adapted for use 
in the present problem. The notation used for the Bessel 
functions is that of Watson (reference 13 ) and of the Smith- 
sonian Tables (reference l 1 *). 


Series for camber effect .— The discussion 
ference effects associated with airfoil camber 
series m 


W 


1 


-y 


e Xa * J, (a b d) 

\ 8 r ) J 1 "(\ 8 r ) 


of the inter— 
involves the 


(A-l) 


convergent for negative values of x. The summation with 
respect to s extends over all the positive roots of the 
eauat ion 

Jj'Us*) - 0 (*2) 


Letting J B = \ 8 r and K= — x, the series may be written 


W .V J .l. ( W r 2 ,-J.K/r- 

1 3 ‘ 


(A3) 


where the summation is taken over all the positive roots of 
the equation 


J i ( J a ) “ 0 


( A 4 ) 


Now, consider the function 
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TT Ji(v) Y/wq/r) - J/wq/r) Y x * (w) -wk/r 
2 j/(w) _ _ 6 


(A 5 ) 


where the quantity Yj, la a Bessel function of the second 
kind of order unity* This function has a simple pole at each 
of the points w = ±J B and Is one-valued and analytic at all 
other points In the ' complex w-plane. Its residue at the point 
j B can he shown to he 


Ja J i(Js) Ji" <J a ) 


-J a K/r 


which is Identical with the general term of the series (A 3 ). 

By the therorem of residues, the integral of the function 
(A5) token counterclockwise around a contour inclosing the 
portion of the complex plane to the right of the imaginary 
axis is then equal to 2-rTiWj.. The integral along a large 
semicircle on tho right of the imaginary axis tends to zero 
when the radiuB of tho semicircle tends to infinity through 
values such that the semicircle avoids the poles of the inte- 
grand. It is thus necessary to retain only the integral along 
the imaginary axis. The contour must, however, have an inden- 
tation to the right of the origin, since the integrand has a 
pole there with residue (r a+qa)/rq. If the radiuB of the in- 
dentation is mado to approach zoro, Wj may finally he written 



IProm the known relations for the modified Bessel func- 
tions, it is readily shown that 


J x (*it) = All/t) 

y (Ait) = -I_(t) A 2 iK (t) 

1 1 rr A 

J 1 '(Ait) = 1/(0 

Ti'^iO - ±il x f (t) + S Z/ (t) 

TT 
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where I x and K x are modified Bessel functions of the 

first and second kind of order unity. By writing the inte- 
gral in equation ( A6 ) in two parts, one along the positive 
and one along the negative imaginary axis, and replacing w 
in these integrals "by -fit and —it, respectively, then 

* ocorc eo 


W, = - 


r 2 + T| s 
3rn 

on 


+ 1 f hhl *ll**h) - ytl/r) 8li (tK/r) dj 

1 1 \ t J 

0 

or co 

W , = — — + — ^ K.(tT)/v) sin(tit/r) dt 


2rp 


1 P K'(t) . . 

- / — — - Iil.tn/r) s in( t k/t ) dt 

n -J 1 1 ( t ) 


( A8) 


The value of the first integral in this equation is given 
by Watson in reference 15 as 

CO 




AiCtTi/r) s in( t k/t ) dt = 


rn-r. 

2T) v / ri a + 


(A9) 


The second integral can be evaluated by expanding the 
product IiCtn/r) sir.(tK/r) in ascending powers of t and 

integrating term by term. The series expansion for the prod- 
uct is 


I ^ tq /r ) r. in( t */r ) 

Y v (-1) P t^ k+p+l) ^sk+i K ap+1 
L L> k! (k+l)i ( 2p+ 1 ) I 2 ak+1 r ^k+P +1 ) 

k = 0 p= 0 


and the t erm— by— t ern integration givos 
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k! ( k+1 ) ! (Sp+1)! 2» k+1 r a t k+p+i ) 

k=0 p=0 

The coefficient ^ k+p+i ) = 2 f 18 siven *7 
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2f " tt J I ' 
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(t) 


which may be written after integration by parts 


co 


2f 


( 2f +1 


i r t af ~ a 

+ 1)tt J [Ix 


af " a (1+t. 3 ) 


(t)] 


dt 


(All) 


This integral is a constant for any given value of f. 

Reverting to the original variable x, the expansion 
for 'a'x may finally be written 


- - 


r s + ri 1 


rx 


2rri 2Tv/Ti a +x a 


°° “ (~D P I*"., % x8p+ i 


. V V * ^ Z 8(to-P+D . (Lia \ 

i kl (k+l)l (2p+l)l 2 ak+1 r a ( k+p+ 1 5 

JC=0 p=o 

This agrees with the result given without derivation by Tani 
and Talma (reference 16). 

t 

For purposes of Computation the coefficient ^ Z f ia 


written 


2f 


( 2f +1 


^^(f-l) + P*2f} (A13) 
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where 


_ne quantity 



& f 2f 


CD 

i r iif. fit 

"J Cii'( t )] s 

o 


( A14 ) 


can then he expressed in a form suitable 


for computation by means of a method devised by Watson for 
an analagous integral (reference 15). 


As the first step, the function 


t 

[Ii'(t)] 3 cosUt/b) 


( A 15 ) 


is written as a sum of partial fractions, b being a posi- 
tive constant which will be fixed later. This can be accom- 
plished by considering the integral 


w ar d w 

J ( w— t ) [Ij w)] s coa(iTw/b) 


around the circle ! w| = H. in the complex plane. The inte- 
grand of the integral (A16) hns poles at the points 


w = t , w = ±ij B , w = ±( n + l/2) b 

where j B is a positive zero of J l , (w); and s = 1, 2, 3, 

. . . , n = 0 , 1 , 2 , 3 , ... . The residue at the simple polo 
at w = t is the function ( A15 ) . The residues at the 
simple poles at w = ±( n + l/2)b are 

_ C-i) n (n + 1/2) af b sf + 1 

tt( nb + 1/2 b?t; [I 2 ‘(nb + l/2 b )] 3 


The poles at w = ±ij a 
there are 


are second order polos; the residues 
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(-l) f J B af+3 

00 Bh(nJ a /l>) 

2f — (nj a /b) tanh( ttJ B /" b) + 


_ j a 

I (j B ^it) a 

a ~J> 1 
l-Js 8 


5 0 ±it 


How, the Integral (A16) taken around the circle |wl = H 
tends to zero when E tends to Infinity In such a manner 
that the circle never passes through a pole of the inte- 
grand. It follows from the theorem of residues that the 
sum of the residues of the integrand at all its poles Is 
zero; thus 


00 


[Ij/Ct)] cosCnt/h) 
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/ i \ H \ 3 f 1 St + 3 
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J B a -t 


( Js®+t a ) s 

(A17) 


By multiplying this equation by cos(nt/b) and inte- 
grating from — » to + 00 , it can be shown with the aid of 
certain integral relations given by Watson (reference 15, 
p. 36) that 


r t sf dt 

J Ci/Ct)] 8 
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and therefore 
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( A19 ) 


i'ho first series in this actuation converges rapidly when b 
is large, the second when b is email. A reasonable com- 
promise for purpoees of calculation is to take b = 1. 


Equation (A19) with b = 1 has been used together with 
equation (AlS) to determine the first four values of the cocf 
ficient * 2 f * flnal results are 


H' a « — C .999 

p' 4 = -1.627 ^ 

n'e = -9.78 ( 

V. \ - - 120.8 ] 


( A20) 


Comparable values of u and p to the same number 

of significant figures are given without derivation in refer- 
ence 18. The value of p ' in this latter reference agrees 
with that of the present paper, but p* 4 differs by one in 
the third decimal place. The value given in (A20) has been 
carefully checked for several values of the parameter b and 
appears to be correct. Values of p"^ and p ' a apparently 
have not previously been computed. 


Series for thickness effect .— The series 
in the discussion of the interference effects 
airfoil thickness is 


W 


e^a* j x ( X b ti ) 

L J»U.r) 

8 = 1 0 8 


which appears 
associated with 


(A21) 
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convergent for negative values of x. The summation with re- 
spect to s ext-end-H J3YGT — all the p-o-a It 1 v e roots of the equa- 
tion w . 

J x (X 8 r) - 0 ( AS2 ) 


Letting J B 
written 


X B r and 


— x ae before, the eerles may be 


V e -J„«/r 

J 2u.) 


/_ 

8 = 1 


(A23) 


where the summation le taken over the positive roots of 


J X (J 8 ) - 0 


( A24 ) 


The function 

n J 1 (w) Y 1 (wq/r) - J i( wn /r ) I^(w) ^_ w */ r 


( A25) 


has a simple pole at each of the points w = ±J a . Its residue 
at eaoh of there polntH can be shorn to be identical with the 
general term of the series (A" r J, Unlike the functinn in the 
previous series, this function is regular at the origin. In- 
tegration around the portion of the complex plane to the right 
o? the imaginary aiis then gives 

, - J_ f 1 n J i(»> T 
a " Sni J 2 

—co i 


Jj/ w) 


By applying the first two of equations (A7) and combining 
the integrals along the two halves of the imaginary axle as 
before, the series becomes 


c* 


= — ■ J' tK^tq/r) cos(tK/r) dt 

QB 

~ ^ S I ~K ~ ) co®(*t/r)dt 

n ” 1' 


(A27) 


The firBt integral can be 
relation (A9) with respect 


evaluated by differentiating 
K j This operation gives 


I 
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/' n r 

j ts^tn/r) eosCtn/r) dt = ^(^a+^a) a/a ( i28 > 

o 

The second Integral can 1)6 evaluated as before by ex- 
panding the preduct tI x (tTj/r) cos(tK/r) in ascending powers 
of t and integrating term by term. The series expanBien 
for the product Is 


tl x (tri/r) oos( tn/r) 


on co 

-II 


(-l) P t a(k + p+l) K 3P 

kl (k+l)I (3p)l 2 ak+1 r ak+ ®P + 1 

k=0 p=0 

and the term— by— term Integration givee 

oo 

ces(tK/r) dt 

. „ f f (-D p i ak+1 

k! ( Jc+ 3.)l(2p)I 2 ak+1 r 2k+ap+i 

where the coefficient ^atk+p+l) = ^af * 8 £* ven 


CD 

/ 


(A29 ) 


1 “ mo 

“•* = - J TJT) * 4t ' 


( 2f +1 


1 f t 3f dt 

+ 1 )t t J I i( t ) 


(A30) 


Reverting to the original variable x, the expansion 
fer W a may finally be written 


w _ r 3 n 

3 " 2(p 3 +x 3 ) 3 / 2 


OO 09 


I I kF 

k=0 p=0 


(-l) P l^a(k+p+l) n^k+i x ap 
(k+1 ) ! (2p)l 2 ak+I r ak+ *P + 1 


( A31 ) 


The integral (A30) has been investigated by Watson (ref- 
erence 15). Its value for any given r can be computed from 
the series 
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4? 


— - * (3f+l) 


sf 


_ _i_ r t af dt ■ 

" J W 

„ i. V ( mu- 1 /a i> ) ,f 

” „4b !;<”> + Vs»> 


+ (-d^ 1 S 


-( 2 f+l) 


v *- 1 

.-i J * S(J * ) 

e -n3 g /Ti 

cosh( nj a /b; 


(nJa/fr) „ 
coBh a (rrj B /b) 


(A3 2) 


where b i 8 an arbitrary positive constant. This equation 
has "been used with b = 1 to determine the first feur values 
of tigf The final results are 


M- e = 0,797 " 
|x 4 = 1.300 
U 6 = 7.46 
H 8 = 96.3 


( A33 ) 


The first two of these values agree ts the three decimal 
pltxce 3 with the two numerical values computed by Watson. The 
remaining two values have not previously been computed. 


APPENDIX B 


CONSTANCY OP LIFT OVEH AN AIREOTL SPANNING 
A CLOSED-THROAT TUNNEL 


Consider an infinitesimally thin untwisted airfoil of 
constant chord spanning a closed— throat wind tunnel of- arbi- 
trary section. Such an arrangement is shown lh figure 10, 
which is a section of the tunnel as seen from downstream. It 
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1b assumed that the flow In the tunnel is nonviscous and that 
the airfoil therefore has no drag. 

Suppose for the time being that the lift varies in some 
manner across the span of the airfoil. Any such variation 
will be accompanied by a system of vortices trailing from the 
airfoil and extending infinitely far downstream. If the usu- 
al assumption is made that the trailing vortices are parallel 
to the axis of the tunnel, the flow pattern in a plane normal 
to the axis at infinity downstream must be of the nature shown 
in figure 10. The flow pattern, in general, consists of a 
number of separate sections within each of which the flow has 
a closed, circulatory character. The line AB, which repre- 
sents the projection of the airfoil, extends across every 
such section, and each of the sections contains the filaments 
of a portion of the system of trailing vortices. The exact 
character of the flow pattern in any particular case depends 
upon the spanwlse variation in lift and upon the cross- 
sectional shape of the tunnel. 

How, consider the flow around a streamline within any 
one of the separate sections of the flow pattern — say the 
streamline OB in the section at the left-hand side of the 
tunnel in figure 10. ThlB streamline, like all the stream- 
lines, intersects the projection AB of the airfoil in two 
points, denoted as C and B in the figure. The fact that - 
in the presence of the tunnel walls each streamline must in- 
tersect AB in two points is essential to the discussion. 

If it is supposed for purposes of discussion that the direc- 
tion of flow is clockwise as indicated, the vertical component 
of velocity at C is upward while the corresponding component 
at B is downward. This direction of flow corresponds to a 
net circulation in the clockwise direction for all the trail- 
ing vortex filaments enclosed within the Btreamllne. 

At the position of the airfoil the pattern of transverse 
velocities induced by the trailing vortices is geometrically 
similar to the pattern at infinity downstream, only the magni- 
tude of the velocities being different. Hence, at points on 
the airfoil directly ahead of point C, the vertical veloc- 
ity induced by the trailing vortices is upward. At points 
directly ahead of point B, the velocity is downward. Thus, 
since the airfoil is untwisted, the airfoil section corre- 
sponding to C operates at a larger effective angle of at— 
taok than does the section corresponding to B, If the air- 
foil is of oonstant chord as assumed, this means that the 
lift at section C must be greater than the lift at section B. 
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As has been pointed out, however, -the trailing vortices 
discharged between sections C and D must have a net cir- 
culation in the olockw-lse direction in figure 10. ' This means 
that the circulation of the spanwlse hound vorticds at sec- 
tion D must he greater than at section C. Since the di- 
rection of stream flow was taken to he toward the observer, 
this in turn means -that the lift at section C must he less 
than that at section '0, which is in direct contradiction 
to the previous result. The original supposition that the 
lift varies across the span thuB leads to two mutually con- 
tradictory conclusions and is therefore invalid. It follows 
that the spanwlse distribution of lift is uniform across an 
untwisted, constant— chord airfoil spanning any closed— throat 
wind tunnel , irrespective of the cross-sectional shape of the 
tunnel . 

As mentioned at the outset, this result depends upon the 
assumption that the airfoil is infinitesimally thin and has 
no drag. It will not he Btrictly true if the Increase in ef- 
fective stream velocity caused by the Interference between 
the walls and the airfoil thickness and wake is not uniform 
across the span. The result also neglectB any effect that 
the boundary layer along the walls of the tunnel may have 
upon the lift distribution. That these appr eximat ions are 
not serious, at least in the case of the circular tunnel, is 
indicated by the experimental results of figures 6 and 7. 

The foregoing reasoning is, of course, Inapplicable for 
an airfoil which does not span the tunnel or for a finite- 
span airfoil in free air. In these instances, the projection 
of the airfoil does not extend across all of the sections in- 
ti? which the transverse flow pattern is divided, and the 
streamlines of this pattern need not intersect the projection 
of the airfoil in two points. Under these conditions a type 
of varying lift distribution can be found whioh does not lead 
to a logical inconsistency. 


AFPBNDIX C 

LIST OB’ IMPORTANT SYMBOLS 

c alrfo.il chord 

d diameter ef circular tunnel 
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radius of circular tunnel 
height of rectangular tunnel 

height of rectangular tunnel equivalent to circular 
tunnel with regard to camber effect 

height of rectangular tunnel equivalent to circular 
tunnel with regard to thickness effect 


1/4 J ; chord— height factor with regard to wake 
effect in rectangular tunnel 
0.331 (V) ; chord— diameter factor with regard to wake 
effect in circular tunnol 


rr a f V'N 
48 \h / 


a 


; chord— height factor with regard to camber 

and thickness effect in rectangular tunnel 
a 


0.289 


(})'• 


chord— diameter factor with regard to 


camber effect in circular tunnel 

f c y 

0.339 ^— j ; chord— d iamet er factor with regard to 

thickness effect in circular tunnel 

factor depending upon shape of base profile (see equa- 
tion (44) and table I) 

angle of attack 

section lift coefficient 

section quart er— chord— moment coefficient 


section drag coefficient 
stream velocity 
dynamic pressure 
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M Mach, number 

E Reynold's number 


x.y,* 
x,u) ,e 
r« 

dr> /a! 

n.t 

* 01^0 

£ 


<F 


v x* . v E 
67' 


*P 
t e 




rectangular apace coordinates 

cylindrical space coordinates (see equations (l)) 

circulation of single line vortex in tunnel 

circulation per unit chord length 

7 and a coordinates of elementary vortex 

radial and angular coordinates of elementary vortex 

chordviae coordinate of elementary vortex; also 

variable of integration in equations (2) and (3) 

velocity potential 

x and s components of induced velocity 

increase in axial velocity at position of airfoil 
in tunnel 

doublet strength 

projected thickness of airfoil 

effective thickness of airfoil 

series of terms involving Bessel functions (Bee 
equations (Al) and (A2l) of appendix A) 


Jm«*m Bessel functions of first and second kind of order 
m (Watson's notation) 

modified Bessel functions of first and second- kind 
of order m (Watson's notation) 

A b variable of summation defined by the roots of the 

equation Jm^Xgr) = 0 

Je A g r ; root of the equation Jm*(Js) “ 0 
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** nf numerical coefficients (see equat ions ( 35 ) and 

(13)) 

k,n, p variables of summation 

t v variables of integration 

K alternalo variable defined as equal to — x 

Superscripts 

(') when pertaining to fluid properties, denotes values 

in the undisturbed stream in the tunnel; when per- 
taining to airfoil character ics , denotes values in 
tunnel, coefficients being referred to dynamic 
pressure q 1 

( # ) denotes firBi derivative of Bessel function with re- 

spect to its argument 

( ' ' ) denotes second derivative of Bessel function with 
respect to its argument 
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TABLE I.- VALUES OP A 



OR VARIOUS BASE PROPILES 


EACi low-drag sections 


H 

<7 

kr 

t 

B 

• t<\ 

B 

S’ 

JP 

<7 

VI) 

a 

3 

VO 

H 

9 

lf\ 

VO 

66-0 xx 

H 

«? 

S' 

0.119 

O.107 

0.106 

0.101 

0.106 

0.109 

0.117 

0.119 

.185 

.155 

.157 

.148 

.158 

.163 

.174 

.178 

.256 

.217 

.214 

.201 

.215 

.221 

.236 

.243 

.330 

.281 

.278 

.257 

.276 

.283 

.301 

.311 

M 

.3US 

.344 

.316 

.338 

.348 

.368 

•379 

.439 

.417 

.415 

.378 

.402 

.414 

.436 

.450 

.602 

.512 

■513 

.460 

.487 

.501 

.526 

.543 








mm 








B 
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gj 
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TABLE II.- COHEHESSIBILIQJT FACTORS FOR CORRECTIO'J EQUATIONS 
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Figure 3.- Infinitesimally Thin Airfoil 5panning 
Closed-Throat Circular Tunnel. 
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Figure 10.- Assumed Flow Pattern m Plane Normal 
to Tunnel Axis at Infinity Downstream 
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